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1 Introduction 



Discrete orthogonal polynomials with respect to uniform lattices have at- 
tracted the interest of researchers from many points of view [25J . A first ap- 
proach comes from the discretization of hypergeometric second order linear 
differential equations and thus the classical discrete orthogonal polynomials 
Charlier, Krawtchouk, Meixner and Hahn appear in a natural way. As a 
consequence of the symmetrization problem for the above second order dif- 
ference equations, you can deduce the (discrete) measure with respect to such 
polynomials are orthogonal. This yields the so-called Pearson equation that 
the measure satisfies. 

In the last twenty years, new families of discrete orthogonal polynomials 
have been considered in the literature taking into account the so-called canon- 
ical spectral transformations of the orthogonality measure. In particular, 
when you add mass points to the discrete measure (Uvarov transformation) 
the sequences of orthogonal polynomials with respect to the new measure 
have been studied extensively (see [12], [3J, [2J, among others). When you 
multiply the discrete measure by a polynomial (Christoffel transformation), 
some results are known [29] . 

From a structural point of view, some effort has been done in order to 
translate to the discrete case the well stated theory of semiclassical orthog- 
onal polynomials (see p3j). In particular, characterizations of such polyno- 
mials in terms of structure relations of the first and second kind, as well as 
discrete holonomic equations (second order linear difference equations with 
polynomial coefficients of fixed degree and where the degree of the polynomial 
appears as a parameter) have been done [21]. Linear spectral perturbations 
of semiclassical linear functionals have been studied in the Uvarov case |17j . 

On the other hand, we must point out that the linear canonical spectral 
transformations (Christoffel, Uvarov, Geronimus) of classical discrete orthog- 
onal polynomials yields discrete semiclassical orthogonal polynomials. But, 
as a first step, it remains open the problem of classification of discrete semi- 
classical linear functional of class one. The symmetric discrete semiclassical 
linear functional of class one have been described in [21]. Notice that the 
classification of D- semiclassical linear functional of class one was done in [7] 
and for class two in [22] . 

The aim of this work is to provide a constructive method of .D^-semiclassical 
orthogonal polynomials based on the Pearson equation that the correspond- 
ing linear functional satisfies. We will focus our attention on the classifi- 



2 



cation of Dx-semiclassical linear functional of class s = 1. In such a way, 
new families of linear functionals appear. Notice that an alternative method 
is based on the Laguerre-Freud equations satisfied by the coefficients of the 
three term recurrence relations associated with these orthogonal polynomi- 
als. Their complexity increases with the class of the linear functional and 
the solution is cumbersome. Basic references concerning this approach are 
H5] as well as [21]. 

The structure of the manuscript is as follows. Section 2 deals with the 
basic definitions and the theoretical background we will need in the sequel. In 
Section 3 we describe the .Di-classical linear functionals, as Di-semiclassical 
of class s = 0. It will be very useful in the sequel taking into account that 
most of the semiclassical linear functionals of class s = 1 are related to them. 
Indeed, in Section 4, the classification of such semiclassical linear functionals 
is given. Some of them are not known in the literature, as far as we know. 
Finally, in Section 5, limit relations in terms of their parameters for semi- 
classical orthogonal polynomials are studied. 

2 Preliminaries and basic background 

Definition 1 Let {fJ> n } n>0 be a sequence of complex numbers and let C be a 
linear complex valued function defined on the linear space P of polynomials 
with complex coefficients by 



Then, £ is called the moment functional determined by the moment sequence 
{/i ra } n>0 and fx n is called the moment of order n. 

Given a moment functional £, the formal Stieltjes function of £ is defined 



For any moment functional £ and any polynomial q(x), we define the moment 
functional q£ by 



(£,x n )=fi. 



by 



oo 




n=0 
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Definition 2 Let C be the linear functional associated with the moment se- 
quence {fi n } n>0 an d 

/J,Q /il • • • /i n 

A n = det . . . 

Then, C is said to be regular/ quasidefinite (resp. positive definite) if A n ^ 
(resp. > 0) for all n G No- 

Definition 3 A sequence of polynomials {P n (x)} n>0 , deg (P n ) = n, is said 
to be an orthogonal polynomial sequence with respect to a regular linear func- 
tional C, if there exists a sequence of nonzero real numbers {Cn}n>o such 
that 

(C, P k P n ) = CrA.n, k, n E N . 

If ( n = 1, then {p n (x)} n>0 is said to be an orthonormal polynomial 
sequence. Notice that if the linear functional is positive definite there exists a 
unique sequence of orthonormal polynomials assuming the leading coefficient 
is a positive real number. 

If the leading coefficient of P n {x) is 1 for every n G Mo, then the sequence 
is said to be a monic orthogonal sequence (MOPS, in short). 

Theorem 4 Let {b n } n>0 and {7n}„>o with 7„ ^ for every n G No be 
arbitrary sequences of complex numbers and let {P n (x)} a sequence of monic 
polynomials defined by the three term recurrence formula (TTRR) 

P n +i(x) = (x- b n ) P n (x) - 7„P n _i(x), (1) 

with P_i = and P = 1. Then, there is a unique linear functional £ such 
that C(l) = 70 and 

(£, P k (x) P n (x)) = 7071 • ■ ■ 7n4,n- 

Proof. See HH Theorem 4.4] ■ 

If the linear functional is positive definite and {Pn(%)} n>0 is the corre- 
sponding orthonormal polynomial sequence, then the above TTRR formula 
becomes 

a n +iPn+i{x) = (x- b n )p n (x) - a„p„_i(x), 
where a n is a real number and a\ = 7„. 
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Definition 5 Let £ be a linear functional and U* : P — > P be a linear oper- 
ator. The linear functional U£ is defined by 

(U£,P) = -{£,U*P} , PeP. 

Example 6 If U is the standard derivative operator D, we have U* = U = 
D. 

Definition 7 A regular linear functional £ is called U-semiclassical if it 
satisfies the Pearson equation U (<fi£) + ipC = or, equivalently, 

(U(<f)C)+ip£,P) = 0, PeP, 

where <p, ip are two polynomials, and <p is monic. The corresponding orthog- 
onal sequence {P n {x)} n>0 is called U-semiclassical. 

In the literature, semiclassical linear functionals with respect to different 
choices of operators have been studied. In particular, ifU = D (the standard 
derivative operator), the theory of P-semiclassical linear functionals has been 
exhaustively studied by P. Maroni and co-workers (see [23] for an excellent 
survey on this topic). 

If U = D u , where 

f(x + u) - f(x) 

D UJ f{x) = , 

cj 

a regular linear functional £ is said to be P^-semiclassical if there exist 
polynomials 4>, ip, where <p is monic and deg0 > 1 such that ((f)£)+ip£ = 0. 
Notice that 

D 1 f(x) = f(x + l)-f(x)=Af(x), 
D_ 1 f(x) = f(x)-f(x-l) = Vf(x), 
are the forward and backward difference operators, respectively, and 

hm PJ(x) = Df(x) = f(x). 

If U — D u , we define U* = D- w . With this definition, we have A* = V and 
when oj — y we recover the identity U* = D = U. 

The concept of class of a P^-semiclassical linear functional plays a central 
role in order to give a constructive theory of such linear functionals. 
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Definition 8 If C is a D u - semiclassical linear functional, then the class s 
of £ is defined by 

s— minmax {deg0 — 2, deg^ — 1} . 

(j>,ip 

among all polynomials 0, if> such that the Pearson equation holds. Notice that 
the class s is always nonnegative. 

For any complex number c, we introduce the linear application 6 C : P — > P 
defined by 

x — c 

We have the following result |23j . 

Theorem 9 The regular linear functional C satisfying the Pearson equation 

D u {4>C) + tpC = 

is of class s if and only if 

11 (|v(c -u) + (e c <p)(c -u)\ + \(c, e c -M + e e <f>)) |) > o, 

cez(4>) 

where Z{4>) denotes the set of zeros of the polynomial 4>(x). 
When there exists c G Z(<j)) such that 

i>{ c -u) + {e c( j>){c -co) = (c, dc-uiti) + e c( j>)} = o, 

the Pearson equation becomes 

[(^)£] + [^(^ + ^ c 0)]£ = O. 

Remark 10 When s = 0, the D^-classical orthogonal polynomials appear 
(see \T$). When to — 1 several characterizations of classical orthogonal poly- 
nomials have been done in flfi]/ . Indeed, we explain with more detail in the 
next section the main characteristics of these polynomials and their corre- 
sponding linear functionals. 

The Di-semiclassical linear functionals have been studied by F. Marcelldn- 
Salto and they are characterized following the same ideas as in the D case. 
Maroni-Mejri deduced the Laguerre- Freud equations for the coefficients of 
the TTRR of D w - semiclassical orthogonal polynomials of class s — 1. In 
the symmetric case, i.e., the moments of odd order vanish, they deduce the 
explicit values of such coefficients and the integral representations of the cor- 
responding linear functionals are given. 
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On the other hand, the Pearson equation yields a difference equation for 
the moments of the linear functional and, as consequence, we get a linear 
difference equation with polynomials coefficients satisfied by the Stieltjes 
function associated with the linear functional. Indeed, 

Theorem 11 If £ is a D^-semiclassical moment functional, then the formal 
Stieltjes function of £ satisfies a non homogeneous first order linear difference 
equation 

<P(z)D ul S c (z) = a(z)S c (z)+b(z), 

where a(z) and b(z) are polynomials depending on <fi and ip, with deg (a) < 
s + 1 and deg(6) < s. 

3 Discrete semiclassical orthogonal polyno- 
mials 

In the sequel, we will consider linear functionals 

oo 

(£, P) = P(x)p (x) , 

x=0 

for some positive weight function p (x) supported on a countable subset of 
the real line. With this choice, the Pearson equation 

(A(0£) + ^£,P) =0, PeP, 

yields 

A(0p)+^p = O. (2) 

We rewrite this equation as 

p(x + l) <f)(x)—ip(x) X(x) 
p(x) = <f>(x + l) = <j)(x + 1)' 

with 

(j) (x) = x (x + ft) (x + (3 2 )---(x + (3 r ), r < s + 2, 

and 

A (x) = c (x + cci) (x + a 2 ) • • • (x + oli) , I < s + 2. 
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Since the Pochhammer symbol {ot) x denned by 



r (a + x) 



satisfies the identity 

( a )x+i 

= x + a, 



Wx 



we obtain 

/ \ Mx ' ' " ( a l)x c 

p(x) - 



(A + 1), • • • (ft + l) x d 
We will denote the orthogonal polynomials associated with p (x) by 

p n r) (x;ai,...,a ; ,/3i,...,/3 r ;c). 
3.1 Discrete classical polynomials 

When s = 0, we solve the Pearson equation (j2j) with deg (0) < 2, deg (ip) = 1 
and three canonical cases appear [25] : 

1. Pn°'°^ (x;c) Charlier polynomials 
We have 

p{x) = — , c>0, xeN , (3) 
(f) (x) = x, if) (x) = x — c, (4) 

and 

X(x) = c. 

The Charlier polynomials have the hypergeometric representation [19] 
C n (x;c)= 2 F (~n, -z; ; ™J (5) 

and the monic Charlier polynomials C n (x; c) become 

C n (x; c) = (-c) n C n (x; c) . (6) 
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2. Pn 1 ' ^ (x; a; c) Meixner polynomials 
We have 

p(x) = (a) x — a > 0, 0<c<l, x E No, (7) 
(x) = x, ^ (x) = (1 — c) x — ca, (8) 

and 

A (x) = c (x + a) . 

The Meixner polynomials have the hypergeometric representation [T9j 
M n (x; a, c) = 2 F X (— n, -x; a; 1 — (9) 

and the monic Meixner polynomials M n (x; a, c) become 

M n (x;a,c) = (a) n (t^j) M n(x;a,c). (10) 

If a = —N, we obtain the Krawtchouk polynomials K n (x; N\ c) where 
p(x) = (-N) t *, c<0, NeN, xe [0,N], (11) 

and 



x X! 



<f>(x)=x, ip (x) = (1 - c) x + cN. (12) 



3. Pn' 1 * 1 (x; a, — N, f3; 1) Hahn polynomials 
We have 

where a, /3 + 1 £ [—AT, 0] , a (/3 + 1) < 0, iV G N, 

0(x) =x(x + /3), ip(x) = ((3-a + N)x + aN, (14) 

and 

A (x) = (x + a) (x - N) . 
The Hahn polynomials have the hypergeometric representation [19] 

Q n (x; a, /3, JV) = 3 F 2 (-n, -x, n - 1 - /3 + a - N; a, —N; 1) . (15) 
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4 Discrete semi-classical polynomials of class 
1 

When s — 1, we solve the Pearson equation ([2]) with deg (0) < 3, deg (tp) < 2, 
and obtain the following canonical cases: 

1. Pn (x; fa, c) Generalized Charlier polynomials 

"^ = Tfih)S- " No ' (16) 

where > — 1 and c > 0. 

2. Pn ^ (x; c) Generalized Meixner polynomials 

where a (/3 + 1) > and c > 0. 

3. Pi 2,0 '' (x; a, —A"; 1) "Generalized Krawtchouk polynomials" 

p{ x ) = {a) x {-N) x - xe [0,iV], 

where a < —N and iVeN. 

4. Pn (x; oti, a 2 ; fa, c) "Generalized Hahn polynomials of type I" 

where «i«2 (/? + 1) > and < c < 1. 

5. P n (x; «i, a 2 , — N, fa, fa, ; 1) "Generalized Hahn polynomials of type 
II" 

P[) {fa + l) x {fa + l) x x\ ' XG [U ' iVJ [U) 

where «i, a 2 ^ [—AT, 0] , fa, fa $ [-N - 1, -1] , a x a 2 {fa + 1) {fa + 1) < 
and A" G N. 

In the next sections, we study these polynomials in detail. 
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4.1 Generalized Charlier polynomials 

We have 

(f) (x) — x (x + f3) , ip (x) = x 2 + (3x — c, 

and 

A (x) = c. 

The first moments are 

/i = c~hp (2^c) r (p + 1) , /n = c^Vj r (/? + 1) , 

where l v (z) is the modified Bessel function of the first kind [261 10.25.2]. 

In [18] Hounkonnou, Hounga and Ronveaux studied the semiclassical 
polynomials associated with the weight function 

p r (x) = j^y, r = 0, 1, . . . . (20) 

When r = 2, they derive the Laguerre- Freud equations for the recurrence 
coefficients and a second-order difference equation. 

In [32] Van Assche and Foupouagnigni also consider (TSUI) with r = 2. 
They simplify the Laguerre-Freud equations obtained in [18] and get 

1 nu n 

y/cU n+1 U n 

with 7„ = c (1 — u^) and (3 n = v n + n. They show that these equations are 
related to the discrete Painleve II equation dPn [31] 

(an + b) x n + c 

•En+l X n —\ - - . 

1 — x n 

They also obtain the asymptotic behavior 

lim 7 n = c, lim v n = 0, 

n— >oo n— >oo 

and conclude that the asymptotic zero distribution is given by the uniform 
distribution on [0,1], as is the case for the usual Charlier polynomials |2Uj . 



u n+1 + U n -i = 

V n = 
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In [50], Smet and Van Assche studied the orthogonal polynomials as- 
sociated with the weight function ( 1161) . They obtained the Laguerre- Freud 
equations 

K +1 - c) (a* - c) = c (6„ -n){b n -n + /3), (21) 

CTi 

b n + 6 n _i = n - 1 - + — , 

for the orthonormal polynomials. They showed that these equations are are 
a limiting case of the discrete Painleve IV equation dPpv [S] 

_ (y n -5n-Ef-A 

vl-B 

5n + E -S/2-C 5n + E — 5/2 + C 
yn + yn - 1 ~ 1 + Dx n + 1 + x n /D ' 

Finally, in Filipuk and Van Assche related the system [21] to the 
(continuous) fifth Painleve equation P v . 

4.2 Generalized Meixner polynomials 

We have 

(f> (x) = x (x + j3) , ip (x) = x 2 + ((3 — c) x — ca 

and 

X(x) — c(x + a) . 

The first moments are 

etc 

fio = M{a,p+ l;c), n 1 = j-^M(a + l,p + 2;c), 

where M(a,b;z) is the confluent hypergeometric function [26l 13.2.2]. 

In [27] Ronveaux considered the semiclassical polynomials associated with 
the weight function 

Pr (x) = Yl {aj) x j^r, r = 1, 2, . . . , 
j=i 

and in [28J he made some conjectures on the asymptotic behavior of the 
recurrence coefficients. 
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In [50], Smet and Van Assche studied the orthogonal polynomials as- 
sociated with the weight function ( 1171) . They obtained the Laguerre- Freud 
equations 



a — 1 . . / a — 1 — (3 \ 

{U n + V n ) {U n+1 + V n ) = —^— v n [v n ~ C) I V n - C - - ^ ) , (22) 



(U n + V n ) [U n+1 + D n _i) = en (« n + C)\u n + C — 

u n a _i \ a i 

for the orthonormal polynomials, with 

= cn — (a — 1) u n , 

On = ^ + a + c- p- l u n . 

c 

They also proved that the system (1221) is a limiting case of the asymmetric 
discrete Painleve IV equation a— dPpy [3T] . 

In [T3] Filipuk and Van Assche showed that the system (1221) can be ob- 
tained from the Backlund transformation of the fifth Painleve equation Py. 
The particular case of (ITT)) when /3 — was considered by Boelen, Filipuk, 
and Van Assche in [9]. 

If we set a = -N, N e N in (ITT]) , we obtain 

(/3 + l^x! 

where we now have /3 > — 1 and c < 0. This case was analyzed by Boelen, 
Filipuk, Smet, Van Assche, and Zhang in [8]. 



4.2.1 Singular limits 

If we let a and (3 — > — 1 in (JIT]) , we have p (x) — » p (#) where p (x) is a 
new weight function satisfying the Pearson equation 

A[(x-l)xp\ + [x-{c+l)]xp = 0. (23) 

Assuming that p(x) satisfies xp(x) = xu(x), for some weight function u(x), 
we get 

A [(x - 1) xu] + [x - (c + 1)] xu = 0. (24) 
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Using the product rule 

A (fg) = fAg + gAf + AfAg (25) 

in ( 12~4"|) . we have 

xu + (x — 1) A (xu) + A (xu) + [x — (c + 1)] xu = 0, 

or 

xA (xu) + [x — (c + 1) + 1] xu = 0. 
Dividing by x, we obtain 

A (xu) + (x — c) u = 0. 

Comparing with (j4]), we see that u(x) is the weight function corresponding 
to the Charlier polynomials ([3]), and therefore (123]) implies that 

p(x) = — + M8(x), (26) 
xl 

where 5 (x) is the Dirac delta function. 

The orthogonal polynomials Pn ^ (x; 0, —1; c) associated with the weight 
function fl2"B"l) were first studied by Chihara in [12]. He showed that they 
satisfy a 3-term recurrence relation ([T]) with 

n £> n , , _ N D n+1 



and 



where 



and 



b n = c—j— + (n + l) , 

n + 1 D n+ i D n 

n 2 Dl 



ln = C 



D 



n 



n + 1 D n _xD n+1 

c n M 
nJ.e c + MK n ^i ' 



K n = Y J - v ^-1 = 0. 



7 

Note that in order that D n is well defined for all n, we need M > — 1, since 
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In [5], Bavinck and Koekoek obtained a difference equation satisfied by 
these polynomials and in [2] Alvarez- No darse, Garcia, and Marcellan found 
the hypergeometric representation 

P&*> (x; 0; -1; c) = (-c) n 3 ^i f-n, -s, 1 + --Y 

Since 

,-►00 (2^ 

we see that 

limP^ 1 ) (ar;0,-l;c) = C n (ar; c) , 

M-s-0 

where C n (x; c) is the monic Charlier polynomial (Q. 

4.3 Generalized Krawtchouk polynomials 

We have 



(ft (x) = x, if) (x) = — c 

and 



x l + [ a- N- - ) r - n.V 



A (x) = c (x + a) (x - N) . 
The first moments are 

/xo = Cn ( — a; — ) , fii = — caiVCjv_i ( — a — 1; — 

V cy v c 

where C n (x; a) is the Charlier polynomial (jSJ). 

To our knowledge, these polynomials have not appeared before in the 
literature. 

4.4 Generalized Hahn polynomials of type I 

We have 

4> (x) = x (x + j3) , 
if) (x) = (1 — c) x 2 + ((3 — coi\ — col<i) x — cai«2, 
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and 

A (x) = c (x + «i) (x + a 2 ) ■ 
The first moments are 

«o = 2-^1 (ai,a 2 ; + l,c) , 

Aii = - - 2-^1 (ai + 1, a 2 + 1; /3 + 2, c) , 

where 2 -Pi («, b; c; z) is the hypergeometric function. 

4.4.1 Singular limits 

If we let ot2 — y 0, /3 — y — 1 and «i = a in ( TT8|) . we have p (x) — )■ p (x) where 
p (x) is a new weight function satisfying the Pearson equation 

A [(x - 1) xp\ + [(1 - c) x - (1 + ca)} xp = 0. (27) 

Assuming that p(x) satisfies xp(x) = xu(x), for some weight function u(x), 
we get 

A [(x - 1) xu] + [(1 - c) x - (1 + ca)] xu = 0. (28) 
Using the product rule fl25|) in (1281 . we have 

xu + (x — 1) A (xu) + A (xu) + [(1 — c) x — (1 + ca)] xu = 0, 

or 

xA (xu) + [(1 — c) x — (1 + ca) + 1] xu = 0. 
Dividing by x, we obtain 

A (xu) + [(1 — c) x — ca] u = 0. 

Comparing with (jSJ), we see that u(x) is the weight function corresponding 
to the Meixner polynomials (j7j), and therefore (|2"5|) implies that 

p(x) = («)^ + M5(x). (29) 

The orthogonal polynomials associated with the weight function ( 129]) were 
first studied by Chihara in [12J. He showed that they satisfy a 3-term recur- 
rence relation ([1]) with 

& _ c (a + n) n g» n + 1 -B n+ i 
c - 1 n + 1 5 n+ i c - 1 S n 
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and 

c n 2 (a + n) B 2 
ln = (c - l) 2 n + 1 B n ^B n+l ' 

where 

g (g) B M 
n (l-c)n!(l- c )- a + Mtf n _ 1 ' 

and 

n j 

Note that in order that B n is well defined for all n, we need M > — 1, since 
^/(l-cp. 

In [10], Richard Askey proposed the problem of finding a second-order 
difference equation satisfied by these polynomials. The problem was solved 
by Bavinck and van Haeringen in [6] , and in [2] Alvarez- No darse, Garcia and 
Marcellan found the hypergeometric representation 

(x; a, 0, -1; c) = (a) n (^T[) ^ 2 ( _n ' ~ X ' 1 + ^ a ' 1 ~ l) ' 
In this case, 

lim P„ (2>1) (x; a, 0, -1; c) = M n (x; a, c) , 

M-»0 

where M n (x; a, c) is the monic Meixner polynomial (1101) . 

If ai = —N, N G N, we can remove the restriction that < c < 1 and 
take any c < 0, with a 2 $ [-N, 0] , (£ [-N - 1, -1] , and a 2 (/3 + 1) > 0. If 
we let «2 — > — (N — 1) and (3 —> —N, we have p (x) —> p(x) where p (x) is a 
new weight function satisfying the Pearson equation 

if) (x) = (1 — c) x 2 + (/3 — c«i — ca 2 ) # — cai« 2 , 

A [x (a: - N) p\ + [(1 - c) x + c (JV - 1)] (x-N)p = Q. (30) 

Assuming that p (x) satisfies (x — iV) p(x) = (x — AT) w(x), for some weight 
function u(x), we get 

A [x (x - N) u] + [(l-c)x + c(N - 1)] (x - N) u = 0. (31) 

Using the product rule ff25|) in fl3Tj) . we have 

xu + (x - JV + 1) A (xw) + [(l-c)x + c(N - 1)] (x - JV) u = 0, 
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or 

(x-N + l)A (xu) + (x - N + 1) (x + Nc - cx) u = 0. 
Dividing by x — N + 1, we obtain 

A (xu) + [(1 - c) x + cN] u = 0. 

Comparing with ffT2~j) . we see that u(x) is the weight function corresponding 
to the Krawtchouk polynomials f lTTj) . and therefore (13"Uj) implies that 

p(x) = (-A0^ + M5 (x — N) . 

4.5 Generalized Hahn polynomials of type II 

We have 

</> (x) = x (x + ft) (x + ft) , 

V> (x) = (ft + ft - at - a 2 + N) x 2 

+ (ftft + aiN - a\a 2 + a 2 N) x + aia 2 N, 

and 

A (x) = (x + «i) (x + a 2 ) . 
The first moments are 

A*o = 3-F2 («i, «2, -iV; ft + 1, ft + 1, 1) , 

^ = - <r ZTVT^ (ai + 1, a 2 + 1, -JV + 1; ft + 2, ft + 2, 1) , 

where 3 F 2 (ai, a 2 , a 3 ; & 1; 6 2 ; z) is the hypergeometric function. 

To our knowledge, these polynomials have not appeared before in the 
literature. 

4.5.1 Singular limits 

If we let a 2 — > 0, ft — > — 1, a\ = a, ft = ft a (fi + 1) < in (|T9|) . we have 
p(x) — )■ p(x) where p(x) is a new weight function satisfying the Pearson 
equation 

A[(x - 1) (x + (3) xp\ + [((3 — 1 - a + N) x + aN — ft xp = 0. (32) 
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Assuming that p{x) satisfies xp(x) = xu(x), for some weight function u(x), 
we get 

A [(x -l)(x + p)xu] + [{p-l-a + N)x + aN- p] xu = 0. (33) 
Using the product rule (125|) in (1331) . we have 

(x + f3)xu + xA [(x + p) xu] + [(p - 1 - a + N) x + aN - p] xu = 0, 

or 

xA [(x + p) xu] + [{P - a + N) x + aN] xu = 0. 
Dividing by x, we obtain 

A [(x + p) xu] + [{p - a + N) x + aN] u = 0. 

Comparing with (I14p . we see that is the weight function corresponding 
to the Hahn polynomials (fIBl) . and therefore (132!) implies that 

?w = ^g*5 +M, w- (34) 

Similarly, if we let a 2 — >■ — (N — 1) , /?2 — > —N, a,\ = a, Pi = P, 
a(P + 1) < in fTl9l) . we have p (sc) — >■ p(x) where p(x) is a new weight 
function satisfying the Pearson equation 

A [x (x + P) (x - N) p\ + [{P - a + N - 1) x + a (N - 1)] (x - N) p = 0. 

(35) 

Assuming that p(x) satisfies (x — N) p(x) — (x — N)u(x), for some weight 
function u(x), we get 

A [x {x + P) {x - N) u) + [{P - a + N - 1) x + a {N - 1)] (x - N) u = 0. 

(36) 

Using the product rule ( 1251) in ( |36j) . we have 

(x + P) xu + (x - N + 1) A [(aj + P) xu] 
+ [(p - a + N - 1) x + a (N - 1)] (x - N) u = 0, 



or 



(x — N + 1) A[(x + P) xu] 
+ (x - N + 1) \{P - a + N) x + aN] u = 0. 



19 



Dividing by x — N + 1, we obtain 

A [(x + P) xu\ + - a + N) x + aiV] u = 0. 



Comparing with (I14p . we see that u(x) is the weight function correspond- 
ing to the Hahn polynomials (ITS]) , and therefore implies that 

p(a;) = ^^i + M5(x-iV). (37) 

The orthogonal polynomials associated with the weight functions ( 134"|) and 
(I37|) were first studied by Alvarez-Nodarse and Marcellan in [I]. 



5 Limit relations between polynomials 

From the two identities [19] 

lim (\a) x = (az) 1 

A— >oo VA/ 



and 

lim — — (Xzf = ( - 

the following limit relations follow: 

1. Generalized Hahn polynomials of type II to generalized Hahn polyno- 
mials of type I 

lim P^ fx; a u a 2 , (3, --, N; l) = P^ (x; a x , a 2 , /3; c) . 

AT-xx> y c J 

2. Generalized Hahn polynomials of type I to generalized Krawtchouk 
polynomials 

lim P^ (x; a, -N, (3; c(3) = P[ 2 >°) (x; a, -N; c) . 

/3— >oo 

3. Generalized Hahn polynomials of type I to generalized Meixner poly- 
nomials 

lim P^ fx; a, a 2 , /?; — J = P^ 1 ' 13 (a; a, /?; c) . 

q 2 ->oo y (^2/ 
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4. Generalized Meixner polynomials to generalized Charlier polynomials 

lim P^) (x;a,ft-) = P<?>» (s;flc). 

5. Generalized Meixner polynomials to Meixner polynomials 

lim P^' l) (x; a, (3; c(3) = M n (x; a; c) . 

6. Generalized Charlier polynomials to Charlier polynomials 

lim (x; (5- cf3) = C n (x; c) . 

We also have the singular limits: 

1. Generalized Meixner polynomials to Charlier-Dirac polynomials 

lim (x; a, (3; c) = C n (x; c) © 6{x), 

a ->■ 
/3->-l 

2. Generalized Hahn polynomials of type I to Hahn polynomials 

lim P^ (x; a, a 2 , (3; c) = Q n (x; a, (3, N) , 
a 2 ->■ —N 
1 

3. Generalized Hahn polynomials of type I to Meixner-Dirac polynomials 

lim Pn^ (x; a, a 2 , (3; c) = M n (x; a; c) © 5(x), 
a 2 ->■ 

4. Generalized Hahn polynomials of type I to Krawtchouk-Dirac polyno- 
mials 

lim P^ l) (x; -N, a 2 , (3; c) = K n (x; -N; c) © 5(x - N), 
a 2 ->• —N + 1 
(3^-N 
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5. Generalized Hahn polynomials of type II to Hahn-Dirac polynomials 

lim P { ^ 2) (x; a, a 2 , —N, f3, f3 2 ; 1) = Q n (x; a, (3, N) © 5(x), 
a 2 ->■ 

6. Generalized Hahn polynomials of type II to Hahn-Dirac polynomials 

lim P T ( t 3 ' 2) (x; a, a 2 , -N, (3, f3 2] 1) = Q n {x; a, (3, N)®5{x - N), 
a 2 ->■ -N + 1 
p 2 -> -N 

where we use the notation "© 5 (x — x ) " to denote the addition of a 
delta function to the measure of orthogonality at the point x . 
We can summarize these results in the following scheme: 



p(3,2) 


->■ 


Hahn© 


5(x) 


Hahn © 5(x — N) 




1 




1 




I 




p(2.1) 
-t n 


->■ 


Meixner ( 


BS(x) 


Krawtchouk © 5(x - N) 


Hahn 


1 


\ 


\ 




I 


1 


p(2,0) 
-t n 


p(l.l) 
-t n 






Charlier © S(x) 






1 


\ 










p(0,l) 
-t n 








Meixner / Krawtchouk 






\ 






1 



Charlier 



6 Concluding remarks 

We have described the discrete semiclassical orthogonal polynomials of class 
s = 1 using the different choices for the polynomials in the canonical Pearson 
equation that the corresponding linear functional satisfies. We have centered 
our attention when the linear functional has a representation in terms of a 
discrete positive measure supported on a countable subset of the real line. 
Some new families of orthogonal polynomials appear as well as some families 
of orthogonal polynomials (generalized Charlier, generalized Krawtchouk, 
and generalized Meixner) which have attracted the interest of researchers 
in the last years taking into account the connection of the coefficients of 
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the three term recurrence relations with discrete and continuous Painleve 
equations. We have also studied limit relations between such families of or- 
thogonal polynomials having in mind an analogue of the Askey tableau for 
classical orthogonal polynomials. It would be very interesting to find the 
equations satisfied by the coefficients of the three term recurrence relations 
for the above new sequences of semiclassical orthogonal polynomials. Fur- 
thermore, an analysis of the class s = 2 will also be welcome in order to get 
a complete classification of such a class as well as to check if new families of 
orthogonal polynomials appear as in the case of the D-semiclassical orthog- 
onal polynomials pointed out in [22] . 



7 Acknowledgements 

The work of the second author (FM) has been supported by Direccion Gen- 
eral de Investigation Cientifica y Tecnica, Ministerio de Economia y Com- 
petitividad of Spain, grant MTM2012-36732-C03-01. 



References 

[1] F. Abdelkarim and P. Maroni. The .D^-classical orthogonal polynomials. 
Results Math., 32(l-2):l-28, 1997. 

[2] R. Alvarez-Nodarse, A. G. Garcia, and F. Marcellan. On the properties 
for modifications of classical orthogonal polynomials of discrete vari- 
ables. In Proceedings of the International Conference on Orthogonality, 
Moment Problems and Continued Fractions (Delft, 1994), J- Comput. 
Appl. Math. (65), pages 3-18, 1995. 

[3] R. Alvarez-Nodarse and F. Marcellan. Difference equation for modifi- 
cations of Meixner polynomials. J. Math. Anal. Appl., 194(1) :250-258, 
1995. 

[4] R. Alvarez-Nodarse and F. Marcellan. The modification of classical 
Hahn polynomials of a discrete variable. Integral Transforms Spec. 
Fund., 3(4):243-262, 1995. 



23 



H. Bavinck and R. Koekoek. On a difference equation for generalizations 
of Charlier polynomials. J. Approx. Theory, 81 (2): 195-206, 1995. 

H. Bavinck and H. van Haeringen. Difference equations for generalized 
Meixner polynomials. J. Math. Anal. Appi, 184(3) :453-463, 1994. 

S. Belmehdi. On semi-classical linear functionals of class s — 1. Classifi- 
cation and integral representations. Indag. Math. (N.S.), 3(3):253-275, 
1992. 

L. Boelen, G. Filipuk, C. Smet, W. Van Assche, and L. Zhang. The 
generalized Krawtchouk polynomials and the fifth Painleve equation. 
larXiv:1204.5W0 vl. 

L. Boelen, G. Filipuk, and W. Van Assche. Recurrence coefficients of 
generalized Meixner polynomials and Painleve equations. J. Phys. A, 
44(3):035202, 19, 2011. 

C. Brezinski, L. Gori, and A. Ronveaux, editors. Orthogonal polynomials 
and their applications, volume 9 of IMACS Annals on Computing and 
Applied Mathematics, Basel, 1991. J. C. Baltzer A.G. 

T. S. Chihara. An introduction to orthogonal polynomials. Gordon and 
Breach Science Publishers, New York, 1978. Mathematics and its Ap- 
plications, Vol. 13. 

T. S. Chihara. Orthogonal polynomials and measures with end point 
masses. Rocky Mountain J. Math., 15(3):705-719, 1985. 

G. Filipuk and W. Van Assche. Recurrence coefficients of a new gen- 
eralization of the Meixner polynomials. SIGMA Symmetry Integrability 
Geom. Methods Appi, 7:Paper 068, 11, 2011. 

G. Filipuk and W. Van Assche. Recurrence coefficients of generalized 
Charlier polynomials and the fifth Painleve equation. Proc. Amer. Math. 
Soc. To appear. 

M. Foupouagnigni, M. N. Hounkonnou, and A. Ronveaux. Laguerre- 
Freud equations for the recurrence coefficients of D u semi-classical or- 
thogonal polynomials of class one. In Proceedings of the VHIth Sympo- 
sium on Orthogonal Polynomials and Their Applications (Seville, 1997), 
J. Comput. Appi. Math. 99: 143-154, 1998. 



24 



[16] A. G. Garcia, F. Marcellan, and L. Salto. A distributional study of 
discrete classical orthogonal polynomials. In Proceedings of the Fourth 
International Symposium on Orthogonal Polynomials and their Appli- 
cations (Evian-Les-Bains, 1992), J. Comput. Appl. Math. 57: 147-162, 
1995. 

[17] E. Godoy, F. Marcellan, L. Salto, and A. Zarzo. Perturbations of discrete 
semiclassical functionals by Dirac masses. Integral Transforms Spec. 
Funct, 5(l-2):19-46, 1997. 

[18] M. N. Hounkonnou, C. Hounga, and A. Ronveaux. Discrete semi- 
classical orthogonal polynomials: generalized Charlier. J. Comput. Appl. 
Math., 114(2):361-366, 2000. 

[19] R. Koekoek, P. A. Lesky, and R. F. Swarttouw. Hypergeometric or- 
thogonal polynomials and their q- analogues. Springer Monographs in 
Mathematics. Springer- Verlag, Berlin, 2010. 

[20] A. B. J. Kuijlaars and W. Van Assche. Extremal polynomials on discrete 
sets. Proc. London Math. Soc. (3), 79(1):191-221, 1999. 

[21] F. Marcellan and L. Salto. Discrete semi-classical orthogonal polynomi- 
als. J. Differ. Equ. Appl, 4(5):463-496, 1998. 

[22] F. Marcellan, M. Sghaier, and M. Zaatra. On semiclassical linear func- 
tionals of class s = 2. Classification and integral representations. J. 
Differ. Equ. Appl, 8:973-1000, 2012. 

[23] P. Maroni. Une theorie algebrique des polynomes orthogonaux. Ap- 
plication aux polynomes orthogonaux semi-classiques. In Orthogonal 
polynomials and their applications (Erice, 1990), volume 9 of IMACS 
Ann. Comput. Appl. Math., pages 95-130. Baltzer, Basel, 1991. 

[24] P. Maroni and M. Mejri. The symmetric .D^-semi-classical orthogonal 
polynomials of class one. Numer. Algorithms, 49(1-4) :251-282, 2008. 

[25] A. F. Nikiforov, S. K. Suslov, and V. B. Uvarov. Classical orthogonal 
polynomials of a discrete variable. Springer Series in Computational 
Physics. Springer- Verlag, Berlin, 1991. 



25 



[26] F. W. J. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark, editors. 
NIST handbook of mathematical functions. U.S. Department of Com- 
merce National Institute of Standards and Technology, Washington, DC, 
2010. 

[27] A. Ronveaux. Discrete semiclassical orthogonal polynomials: general- 
ized Meixner. J. Approx. Theory, 46(4):403-407, 1986. 

[28] A. Ronveaux. Asymptotics for recurrence coefficients in the generalized 
Meixner case. J. Comput. Appl. Math., 133(l-2):695-696, 2001. 

[29] A. Ronveaux and L. Salto. Discrete orthogonal polynomials- 
polynomial modification of a classical functional. J. Differ. Equ. Appl, 
7(3):323-344, 2001. 

[30] C. Smet and W. Van Assche. Orthogonal polynomials on a bi-lattice. 
Constr. Approx. 36: 215-242, 2012. 

[31] W. Van Assche. Discrete Painleve equations for recurrence coefficients of 
orthogonal polynomials. In Difference equations, special functions and 
orthogonal polynomials, pages 687-725. World Sci. Publ., Hackensack, 
NJ, 2007. 

[32] W. Van Assche and M. Foupouagnigni. Analysis of non-linear recurrence 
relations for the recurrence coefficients of generalized Charlier polyno- 
mials. J. Nonlinear Math. Phys., 10(suppl. 2):231-237, 2003. 



26 



